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Abstract. The Feshbach-Villars equations, like the Klein-Gordon equation, are 
relativistic quantum mechanical equations for spin-0 particles. We write the Feshbach- 
Villars equations into an integral equation form and solve them by applying the 
Coulomb-Sturmian potential separable expansion method. We consider bound-state 
problems in a Coulomb plus short range potential. The corresponding Feshbach-Villars 
Coulomb Green’s operator is represented by a matrix continued fraction. 
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I. Introduction 

It is well known that the basic relativistic equation for a spin-0 particle is the Klein- 
Gordon equation. Probably, it had been formulated hrst by Schrodinger, but it gave 
wrong results for the hydrogen atom, so he abandoned it in favor of the non-relativistic 
formalism. The Klein-Gordon equation was suggested in 1926 by Klein [1] and Gordon 
[2] as well as by several other people [31 H [5]. It can be derived from the relativistic 
energy-momentum relation using the correspondence principle. It has served as a basic 
relativistic equation for spin-0 particles in relativistic quantum mechanics and held 
theory. 

The Klein-Gordon equation, however, contradicts some of the postulates of 
quantum mechanics. In the quantum mechanics it is postulated that the system is 
completely determined by the wave function and the time evolution of the wave function 
is determined by the time-dependent Schrodinger equation iniE|. The Klein-Gordon 
equation contains a second derivative in time, therefore to determine the system the 
wave function in not enough, we need its time derivative as well. So, obviously, the 
Klein-Gordon equation is not a genuine quantum mechanical equation. 

In 1958 Feshbach and Villars rewrote the Klein-Gordon equation into a Hamiltonian 
form [8]. They split the Klein-Gordon wave function into two components and for 
the components vector they arrived at a Schrodinger-like equation with hrst order in 
time derivative. Although the Feshbach-Villars formalism appear in some advanced 
quantum mechanics books p cni ni ca [i3i El [15], and they were utilized in gaining 
deeper insight into relativistic physics of Klein paradox pair production, UniElEl 
flp. [201 [211 ISa [23], in exotic atoms [231 [25l used in theoretical consideraions 
[27] [28] EHl El], study relativistic scattering [32l [32] and optics [33] or demosntrate PT 
symmetry [331ES ESI EZ] , they were hardly used as a computational tool. The equations 
look like ordinary coupled differential equations, but the components are coupled by the 
kinetic energy operator, which makes them very hard to solve. 

The aim of this work is to develop a solution method to relativistic quantum me¬ 
chanical problems by solving the spin-0 Feshbach-Villars (FVO) equations. We adapt a 
solution method which worked extremely well in the non-relativistic case. First, in Sec. 

II, we outline the FVO formalism. Then in Sec. Ill we recapitulate the solution method 
for the non-relativistic problems with Goulomb-like potentials. In Sec. IV we present 
our approach to the FVO equation. As an example, we choose an attractive Goulomb 
plus short-range potential problem. This also provides us a convenient test case since 
the Klein-Gordon hydrogen levels are well known from standard text books. Finally we 
summarize our hndings and draw some conclusions. 
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2. Feshbach-Villars equations for spin-0 particles 


The Klein-Gordon equation with a scalar potential V is given by 

2 




In the FVO formalism the wave function is split into two components 

T = 0 + x, (^ih^-V^^ = mc^{(j)-x), 

and the components satisfy the coupled equations 
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and the Hamiltonian 
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equation 


p2 

-11 


1 0 
0 -1 


me -f 


1 0 
0 1 


V 


( 6 ) 


and (01) into a form analogous to the time-dependent Schrodinger 


d 

= Hpvol 


or, for stationary states we may have 

Hpvolf^) = Elf:) 


(7) 


( 8 ) 


This Hamiltonian looks like a usual coupled-channel Hamiltonian. However, there 
the channels are coupled by some short-range potential, while here the coupling is due 
the kinetic energy operator. As we are going to see below, the matrix elements of the 
kinetic energy operator on a discrete basis representation behaves like n for large n. 
This is a long-range coupling, it cannot be truncated to a hnite basis representation. 


3. Solution method for non-relativistic problems 

We consider a Coulomb-like potential 


Vl^'fP + vP, 


( 9 ) 
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where I is the angular momentum, = Ze^/r is the Coulomb potential and is 
some short-range potential. The bound states are the solutions of the homogeneous 
Lippmann-Schwinger equation 

\,p)=g^{E)v[^\fj), ( 10 ) 


where 

{E) = {E — /2m — v^)~^ (11) 

is the the Coulomb Green’s operator. 

If we approximate the short-range potential by a hnite-rank operator, we can solve 
the equations without any further approximation, provided we can calculate the matrix 
elements of the Green’s operator. In this respect the use of the Coulomb-Sturmian (CS) 
functions 

(rlii) = ( („_,_2i + l)! j e->-{2brrnp\2br), (12) 

where L is the Laguerre polynomial and 6 is a parameter, is particularly advantageous. 
The CS functions form a basis. With {r|fi) = (r|n)/r we have: 


{n\n') = 


(13) 


and 
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The simple form of CS basis allow an exact and analytic calculations of the matrix 
elements: 


{n\n') 


{n-\-l l)/b 
-^/n'{n' + 2l + l)/(26) 
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if n = n', 
if n' = n -|- 1, 
otherwise. 
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and 


(n 



{n-\-1 -\-1) b/ {2m) 
\/n'{n' 21 1) 6/(4m) 

0 


if n = n', 
if n' = n -f 1, 
otherwise. 


(16) 


Thus the operator J = E—p^/2m — v^ is an inhnite tridiagonal matrix on the CS basis. 
An important result of Refs. [381 139] is that for inhnite tridiagonal matrices the N x N 
part of the inverse, or the Green’s matrix, is given by 


9^=U 


diNdjNJN,N+lCN+lJN+l,N) 


-1 


(17) 


where g^ and J are N x N matrices and G is a continued fraction. So, basically, g^ 
is almost the inverse of J, only the right-down corner of J is modihed by a continued 
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fraction. The continued fraction is constructed from the higher-index elements of J and 
is defined by the recursion relation 

Cn+1 = {Jn+1,N+1 — ■Jn+1,N+2CN+2Jn+2,N+i) (18) 

Over many years various finite-rank expansion methods for approximating have 
been used. In a recent work we proposed a simple, straightforward, yet very efficient 
approximation scheme m. It amounts of representing the short-range potential in a 
larger basis, invert the potential matrix, truncate to a smaller basis, and then invert 
it back. This way we achieve a low-rank representation of the potential operator that 
contains the relevant information from the larger basis. Even a low-rank representation 
gives very good figures, while higher-rank representations give extremely accurate 
results. 

This method of solving problems of non-relativistic quantum mechanics through 
integral equations has been extended further for solving more challenging problems like 
the three-body Faddeev integral equations with Coulomb potentials [HI 02] . 


4. Relativistic spin-0 particles in Conlomb-like potentials 


We consider a relativistic spin-0 particle in a Coulomb-like potential. 
Hamiltonian is given by 

Hfvo = HpYQ -f VpYQ 


with 


and 



Its FVO 

(19) 

( 20 ) 


( 21 ) 


Analogous to the non-relativistic case, we can cast the eigenvalue problem of the two- 
component wave function into a Lippmann-Schwinger form 


IV') ~ ffFVo(-^)'^FVol'^) ’ 


( 22 ) 


where 

9™(S) = [J!-ff?v„l-‘. (23) 

Since Vcyi) is just a diagonal 2x2 matrix with in the diagonal, we can 

(s) 

approximate n) ^ as before in the non-relativistic case. 

In the CS basis the constant, the Coulomb potential and the kinetic energy 
operators are either diagonal or tridiagonal infinite matrices. So, for 


d^FVO — 



1 A ^ 

— 1 J 2m 




mc^ — 




,c 




(24) 
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the infinite tridiagonal structure is superimposed by some 2x2 matrix structures. In 
other words, Jpvo is infinite block-tridiagonal with 2x2 matrix blocks. The same 
arguments which resulted in Eq. (ITTll is applicable here and leads to a similar expression, 
only Jjj’s should be understood as 2 x 2 matrices and the continued fraction should be 
replaced by a continued fraction of 2 x 2 matrices, i.e. the continued fraction becomes 
a matrix continued fraction. 

In CS basis representation fl2^ turns into a homogeneous algebraic equation 

(25) 

which is solvable if the determinant is zero 

( 26 ) 


5. Numerical illustrations 

In this work we use units such that h = 1, m = 1, = 1, and c = 137.03602. We 

consider first the hydrogen atom with Z = —1 and / = 0. The table below show 
the exact non-relativistic and relativistic Klein-Gordon hydrogen energies, as well as 
the numerical results of the FVO equations. We found that the numerical zeros of the 
determinant in fl26l) coincide with the exact Klein-Gordon hydrogen eigenstates up to 
10 digits, even for higher excited states, irrespective of N. 

Then we add a short-range potential 

= Vo e-"“7r (27) 

with vq = 2 and oq = 2, and investigate its effect on the bound states. The results 
are shown in the last column. We found that the method is accurate enough to give 
account of the fine relativistic effects. We also experienced a rapid convergence in N, 
the number of states used in the expansion. The rate of convergence was similar that 
we observed in the non-relativistic case in Ref. [lO] . 

6. Summary and conclusions 

In this work we considered a relativistic spin-0 particle in a Goulomb plus short-range 
potential. We adopted the Feshbach-Villars formalism. In this formalism the Klein- 
Gordon wave function is split into two components and for the component we get a 
Schrodinger-like equation with a non-hermitian Hamiltonian. We have rewritten the 
equations into a Lippmann-Schwinger form and represented the short-range part of the 
potential in GS basis by a finite matrix. We calculated the GS matrix elements Feshbach- 
Villars Goulomb Green’s operator by a matrix continued fraction. This representation 
of the Green’s operator is exact, thus even a small basis representation gives a faithful 
account of the whole spectra of the Klein-Gordon hydrogen atom. The method gives 
very accurate results for the Goulomb plus short-range case as well. 
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n 

Sch-t-n*^ 

KG+n^ 

FVO+n^ 

FV0-l-n‘"-|-n® 

0 

-0.50000000 

-0.50003329 

-0.50003329 

-0.28203629 

1 

-0.12500000 

-0.12500541 

- 0.12500541 

-0.09299089 

2 

-0.05555556 

-0.05555728 

- 0.05555728 

-0.04546869 

3 

-0.03125000 

-0.03125075 

-0.03125075 

-0.02685280 

4 

-0.02000000 

-0.02000039 

-0.02000039 

-0.01770225 

5 

-0.01388889 

-0.01388912 

-0.01388912 

-0.01254053 

6 

-0.01020408 

-0.01020423 

-0.01020423 

-0.00934633 

7 

-0.00781250 

-0.00781260 

-0.00781260 

-0.00723344 

8 

-0.00617284 

-0.00617291 

-0.00617291 

-0.00576369 

9 

-0.00500000 

-0.00500005 

-0.00500005 

-0.00470027 

10 

-0.00413223 

-0.00413227 

-0.00413227 

-0.00390614 


Table 1. Feshbach-Villars eigenstates of hydrogen atom (FVO+t!*^), hydrogen plus 
short range potential (FVO+u*^ + u®) compared to the exact Klein-Gordon (KG+u‘^) 
and Schrodinger (Sch+u*^) results. 


This is another example of non-hermitian quantum mechanics |13]. The FVO 
Hamiltonian is not hermitian, yet it has real eigenvalues. It is hermitian in a generalized 
sense [iniiii]. Interestingly, Feshbach-Villars equations can also be derived by requiring 
CPT symmetry for the quantum mechanical equation |13] . 

The solution of the Feshbach-Villars equations pose a real challenge to most of 
the solution methods due to the long-range-type kinetic energy in the coupling term. 
However, we can wrap up the infinite block-tridiagonal matrix into the Green’s matrix 
by using matrix continued fractions thus eliminating the main hurdle of working with 
Feshbach-Villars equations. 
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